Muonium emission into vacuum from mesoporous thin films at cryogenic temperatures 
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We report on Muonium (Mu) emission into vacuum following fi + implantation in mesoporous thin 
Si02 films. We obtain a yield of Mu into vacuum of (38±4)% at 250 K temperature and (20±4)% 
at 100 K for 5 keV fi + implantation energy. From the implantation energy dependence of the Mu 
vacuum yield we determine the Mu diffusion constants in these films: fl™ = (1.6 ± 0.1) x 10 -4 
cm 2 /s and D™ K = (4.2±0.5) x 10 -5 cm 2 /s. Describing the diffusion process as quantum mechanical 
tunneling from pore-to-pore, we reproduce the measured temperature dependence ~ T 3//2 of the 
diffusion constant. We extract a potential barrier of (—0.3 ± 0.1) eV which is consistent with our 
computed Mu work-function in Si02 of [—0.3, —0.9] eV. The high Mu vacuum yield even at low 
temperatures represents an important step towards next generation Mu spectroscopy experiments. 



Muonium (Mu), the bound state of a positive muon 
(/i + ) and an electron, is a pure leptonic atom. It is thus 
an ideal object for testing bound-state quantum electro- 
dynamics (QED) free from hadronic uncertainties related 
to the structure of the nucleus [H, . A renewed interest 
in this simple system has been triggered by the recent 
results of the muonic hydrogen experiment 3] . The puz- 
zling proton radius discrepancy observed there could be 
ascribed to problems related either to bound-state QED 
theory, or the Rydberg constant, or the structure of the 
proton or new physics. Spectroscopy of Mu addresses the 
first two items. Furthermore, Mu spectroscopy provides 
precise determination of other fundamental constants like 
the muon mass and the fine structure constant dm. Mu 
can also be used to search for new physics such as lepton 
flavor violation via muonium-antimuonium oscillation [7( . 

For next generation experiments, it is essential to have 
a source of Mu with high vacuum yield down to low tem- 
perature and long term stability. Mu in vacuum is typ- 
ically produced by stopping a low momentum fi + beam 
close to the surface of tungsten foils [8] or silica pow- 
ders . The fraction of Mu which diffuses to the surface 
is emitted into vacuum. Prior to this study, the highest 
measured vacuum yield was (18 ± 2)% per stopped 
obtained in SiC>2 powders at 300 K [9t-llll]. Moreover, 
to our knowledge, Mu emission into vacuum below room 
temperature has never been reported. A Mu source with 
a larger flux can be achieved either by improving the /z + 
beam (smaller phase space, low energy, high intensity) 
as proposed in 12 , [l3| or by improving the /i + —> Mu 
conversion. In this work, we focus on the optimization 
of the latter using SiC>2 porous films (F-samples of [lit ]) 
which we preselected with the ETH Zurich slow positron 
beam. The choice of this material was motivated by the 
fact that Positronium (Ps, the electron-positron bound 
state) and Mu share similar formation mechanisms. Re- 
cently, a yield of Ps into vacuum as high as 40% from 



these porous samples has been measured down to cryo- 
genic temperatures [l4j |. 

For this study, we used the low energy positive muon 
beam (LEM) at PSI delivering approximately 3000 s _1 
li + on target with energies tunable from 1 to 30 keV 
[H El]. The /i + are implanted in the porous film of 
1 /xm thickness, pore size of (5±0.5) nm and density of 
1.1 g/cm 3 . The mean implantation depth is 75(270) nm 
for a /i + implantation energy of 5(19) keV. The Mu for- 
mation mechanism is similar to the one in SiC>2 powders 
0, [l3]- The fi + implanted at keV energy in the SiC>2 
film rapidly thermalize in the bulk (in tens of ps). A 
fraction of them forms Mu in the bulk. Those atoms dif- 
fuse until they are ejected in the pores with almost 100% 
probability [9j. The porous films have a network of in- 
terconnected pores in which Mu can diffuse and loose its 
energy via collisions with the pore walls. If Mu reaches 
the film surface before decaying, it is emitted into vac- 
uum. We define the Mu vacuum yield as the probability 
of Mu emission into vacuum per implanted /x + . If Mu 
suffered a sufficient number of collisions, during its dif- 
fusion to the surface, it becomes thermalized at the film 
temperature. While Ps from similar films is emitted into 
vacuum with an energy above room temperature 141 |18| 
due to quantum mechanical confinement in the pores, 
for Mu one does not expect such a limitation because 
the de Broglie wavelength is about 10 times smaller, i.e., 
of the order of 0.4 nm. 

The LEM is a dedicated facility for /xSR (muon spin 
rotation) measurements. A sketch of the sample reg ion 
and the positron detectors is shown in Fig. [T] (see 15j, |l6[ 
for more details). Before stopping in the sample the /i , 
which are almost 100% transversely polarized, cross a 
10 nm thin carbon foil causing the emission of secondary 
electrons. These electrons detected by a micro-channel 
plate provide the event trigger. Segmented plastic scin- 
tillators surrounding the sample region in a cylindrical 
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FIG. 1. LEM sample chamber. The sample is glued on a sil- 
ver coated copper mount contacted to a cryostat. The sample 
is surrounded by scintillators for positron detection grouped 
in upstream and downstream counters. Each of them is addi- 
tionally segmented in top, bottom, left and right. 



geometry are used to detect the positron from muon de- 
cay. The positron signal provides the stop time of the 
event. The sample resides in a magnetic field transverse 
to the muon spin (see Fig. [T]). Therefore, the muon spin 
undergoes Larmor precession whose frequency depends 
on the local magnetic field and on whether the muon re- 
mains a free fi + or binds with an electron to form Mu. 
Since the positron from muon decay is emitted preferen- 
tially along the muon spin, using a segmented detector 
divided in four sections (top, bottom, left, right), it is 
possible to track the spin precession. The time spectra 
measured in each individual segment follow the exponen- 
tial muon decay distribution, modulated at the Larmor 
frequency. The number of counts N(t) measured in one 
of the positron detectors, e.g., upstream top, is [9|: 

N(t) = N e~*/ T [l + A M {t) + A Uu (t)} + B 

where No is the normalization, r = 2.2 lis is 
the muon lifetime, B the uncorrelated background, 
A^t) = A M e~ A "' cos (ui^t - 4>n) and A Mu (t) = 
^M U e~ AMu * cos (ujMut — 4>Mu) are the precession signals at 
frequencies uj^ for free /i + and ojmu for Mu and phases </> M 



and 

the damping of the precession signal amplitudes A^ and 



!>Mu- The constants A M and Amu take into account 



Amu due to spin relaxation processes [17[ • Since the gyro- 
magnetic factor of Mu in the triplet state (F=l, M=±l) 
is 103 time larger than the gyromagnetic factor of /z + 
(wmu ~ 103 0^+), it is possible to clearly distinguish if 
an implanted fi + remains unbound or forms Mu. 

The initial fraction of Mu formed in the sample per 



implanted /i + is determined with F^ u 



= 1 - A^+JAtot 



where A to t is the total observable asymmetry which has 
been measured using a Suprasil (fused quartz) reference. 
The correctness of this indirect approach relies on the fact 
that fi + is not expected to depolarize in silica (ljj . and 
thus the missing /x + fraction is the one that converted to 
Mu. The measured Mu formation probability (see Fig. ^ 
for porous SiC>2 is F^ u — (60 ± 2)% which is compara- 
ble with the results obtained in silica powders @- For 
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FIG. 2. Fyi u versus temperature for the porous film (8 x 10 6 
events) and Suprasil (2.5 x 10 6 events) for various implanta- 
tion energies obtained from the /iSR amplitudes and Amu- 



Suprasil we obtained — (80 ±4)% in agreement with 
[20j] . In the same plot, we show the initial fraction of Mu 
extracted directly using F^ u — 2^4 Mu /^4 tot [§]. As one 
can see, these values differ from the ones obtained in- 
directly from A^. This is because the direct method is 
sensitive only to the fraction of Mu that does not undergo 
fast relaxation, e.g., due to spin exchange collisions in the 
pores Note that because of the large gyromagnetic 
ratio Mu is much more sensitive to depolarizing sources 
than li + . 

Using the standard /iSR setup allows us to determine 
the probability to form Mu. However with this technique 
we are unable to demonstrate Mu emission into vacuum. 
One possibility would be to use a tracking detector as 
in [10( . We developed a new approach which exploits the 
existing liSR setup. The principle is based on the fact 
that the detection efficiency in the downstream detec- 
tors (see Fig. Q]) is time dependent in case of Mu emis- 
sion into vacuum. Positrons from Mu decaying outside 
of the film have a higher probability to be detected in the 
downstream counters than the ones coming from /i + /Mu 
decays in the sample which are shielded by the copper 
sample support. Therefore, if vacuum emission occurs, 
a deviation from the /i + exponential decay distribution 
is expected. Hereafter, we will refer to this method as 
positron shielding technique (PST). Note that in PST we 
do not consider top, bottom, left and right counters sep- 
arately as in the /zSR setup, but we only distinguish be- 
tween upstream and downstream detectors. In Fig. [3] (a), 
we show the time spectra expected in the downstream 
counters from simulations using Geant4 2l| for 0% (/o) 
and 100% (fwo) Mu yield in vacuum. In Fig. [3] (b), we 
present the measured data for Suprasil (no emission into 
vacuum, thus corresponding to 0%) and for SiC>2 porous 
material where emission into vacuum is expected. In or- 
der to determine the fraction of Mu emission into vacuum 
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FIG. 3. (a) Simulated time distributions in the downstream 
detector for 0% (dashed green) and 100% (solid red), (b) 
Measured time spectra for the porous material (dotted black) 
and the Suprasil sample (solid blue), (c) Data and /fit(i) after 
subtraction of the properly normalized exponential muon de- 
cay distribution. The inset shows the prompt peak and fat (t) 
(without subtraction). 



(-^Mu)' we ^ measured time spectra with 

/ fit (t) = n[(l - FXJ/oW + F v M J 10 o(t)} + n pp f pp (t) 

where n is the normalization and n pp f pp (t) accounts for a 
prompt peak. This prompt peak which occurs in the first 
bins of the time spectra (see inset of Fig. 0(c)) originates 
from /i + decaying in flight before reaching the target 
and from back-scattered /i + [22]. The time distribution 
of this peak f pp (t) is determined experimentally using 
the Suprasil sample. The three free parameters of the fit 
are tz, n< pp and FJ, U . Fits of fat{t) to the experimental 
data which have been taken for various implantation en- 
ergies and film temperatures typically give a reduced \ 2 
of 1.1-1.4 (612 degrees of freedom). In the simulations, 
Mu is assumed to be emitted from the surface of the 
sample with a cos 9 angular distribution 



11, IS! and an 



energy spectrum corresponding to a Maxwell-Boltzmann 
distribution at the target temperature. Fitting the data 
with an isotropic angular distribution or a different tem- 
perature worsens the reduced y 2 by more than 0.2. 
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FIG. 4. (a) Vacuum yield FJju versus temperature deter- 
mined with PST. For comparison, we show the results of the 
direct method, (b) -Fmu versus the implantation energy. The 
curves are fit with the diffusion model described in the text. 



In order to better visualize the comparison between 
simulations and measurements, in Fig.[3](c), we show the 
time spectrum after subtraction of the prompt peak and 
the exponential muon decay distribution. The Suprasil 
data give a constant value as expected due to the absence 
of Mu emission from this sample. On the contrary, for 
the porous film there is a clear signal caused by the in- 
creased positron detection efficiency when Mu is emitted 
into vacuum. The values of F\ An extracted from the fits 
are presented in Fig. 2) We obtain a yield of Mu into vac- 
uum of = (38 ±4)% at 250 K and = (20 ±4)% 
at 100 K for 5 keV implantation energy. The abrupt 



change of F^ An visible between 75 K and 100 K is due to 
thermal absorption of Mu at the pore walls as already re- 
ported for silica powders 23|, |24| . For 20 K Fl Au is com- 



patible with zero. The linear dependence of oc T 
(between 100 K and 250 K) is interesting since from a 
classical diffusion model a V^F-dependence is expected. 
For comparison in Fig. [4] (a), we present also the fraction 
of polarized Mu determined directly from the measure- 
ment of the Mu asymmetry Ayi u with /iSR technique. As 
one can see, these points are systematically lower than 
Fy u obtained with PST. This is because PST, in con- 
trast to the /xSR direct method, is also sensitive to the 
fraction of Mu that depolarizes fast. Nevertheless, both 
methods give consistent results in terms of dependence 
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on the sample temperature and p + implantation energy 
(E). 

In Fig. I] (b), F v Mn versus E at 100 K and 250 K 
is fitted with a one dimensional diffusion model origi- 



nally developed for Ps [25J, |26j. The Mu fraction dif- 
fusing into vacuum is given by F^ u (i?) = F^ lu (E)J(E) 
with J(E) = J Q e~@ x P(x, E)dx, where I is the him 
thickness, j3 = 1/\ZDm^t the inverse of the diffusion 
length and £>Mu the diffusion coefficient. For the ini- 
tial Mu fraction F^ u we used the point at E = from 
Fig. [2j The p + implantation profile P(x, E) was cal- 
culated using the TrimSP simulation validated for p + 
with experimental data (27j . The only fit parameter to 
the data is the Mu diffusion constant £>Mu- The re- 
sulting values determined from the fits (solid lines in 
Fig. I (b)) are D™ K = (1.6 ± 0.1) x 10~ 4 cm 2 /s and 
D™° K = (4.2 ± 0.5) x 10~ 5 cm 2 /s. The good agree- 
ment between fit and data implies that -Dmu does not 
depend on the implantation energy. This means that the 
Mu thermalization time is much shorter than the diffu- 
sion time (cfr. this result with similar measurements in 
Ps, see Fig. 10 of [18]). A further argument that Mu 
quickly thermalizes is given by the worsening of the x 2 
when fitting the data of Fig. [3] (c) with distributions sim- 
ulated at temperatures different from the sample temper- 
ature. Therefore, we can write the diffusion coefficient as 
a function of the mean kinetic energy Em u of thcrmalized 
Mu in the pores as Dmu = A/ (3C) 1 /2i?Mu/'7iMu where 
C is the mean number of collisions that Mu undergoes 
in one pore before reaching the next one, mMu the Mu 
mass and A the mean distance between the pores [lij . 
Assuming an hexagonal close packing of the pores, one 
can estimate the mean separation between them using 
p = p (l - 7rd 3 /A 3 VT8) [Hj]. We obtain A = 5.6 nm 
for a pore diameter of d = 5 nm, a silica bulk density 
of po = 2.2 g/cm 3 and a porous film density of p = 1.1 
g/cm 3 . Using the experimentally determined -Dmu, the 
mean number of collisions in each pore is C = 2100 ± 500 
at 100 K and C = 850 ± 100 at 250 K. These values con- 
firm that Mu thermalization is fast (~ ns) on time scale 
of the diffusion process (~ us). In fact, from the mass 
difference of Mu and SiC>2 [28| , one expects that in order 
to reach thermal energy Mu needs ^500 collisions. 

The obtained Z?Mu values are three orders of mag- 
nitude smaller than expected from a classical diffusion 
model 9]. In order to explain this disagreement, we 
interpret the Mu diffusion process in the porous ma- 
terial as quantum mechanical tunneling from pore to 
pore through a step potential barrier of (0.6±0.2) nm 
width (corresponding to the pore walls thickness) . From 
the mean number of collisions C, which is the inverse 
of the pore-to-pore tunneling probability, we deduce a 
height of the potential barrier of (0.3±0.1) eV. The un- 
certainty is dominated by our poor knowledge of the ma- 
terial structure. With this quantum mechanical model 



we can reproduce the observed dependence of Dm u (T) 
versus the temperature T. Since in our regime the tun- 
neling probability scales approximately linearly with T, 
we obtain that D Mu (T) cx VEy^/C oc Vt/T" 1 cx T 3 / 2 . 
The measured ratio Dj%° K /D$° K » 3.8 ± 0.5 com- 
pares well with the expected value from the T 3 / 2 de- 
pendence of (250K/100K) 3 / 2 R3 4 supporting the va- 
lidity of our model. To check if the value of the po- 
tential barrier height obtained above can be identified 
with the Mu work- function (W), we performed Den- 
sity Functional Theory calculations within GAUSSIAN 
98 [29| on clusters of SiO*2 containing up to eight silicon 
atoms and terminated by oxygen, capped with hydrogen 
atoms. We compute the total energy £^°o 2 +Mu of the 
SiO*2 matrix with a Mu atom and the total energy of 
the SiO*2 fragment alone E^°q 2 . These computations of 

w = ^Si02+Mu - E sTo2 ~ 13 - 6 eV y ield a value between 
-0.3 eV and -0.9 eV. The spread of the interval for W 
originates from the uncertainty to locate the exact posi- 
tion of the interstitial Mu site with respect the Si and O 
atoms. Considering the over-simplification of our model, 
we conclude that our experimental determination of the 
work-function is consistent with the theoretical estima- 
tion. Further experiments using other techniques and 
more precise measurements for Mu and Ps will be useful 
to gain a deeper understanding of this intriguing diffusion 
process in mesoporous films. 

Summarizing, we have found that a sizeable fraction 
of thermalized Muonium is emitted into vacuum from 
mesoporous thin SiO*2 films. At 250 K the yield is more 
than a factor two higher than previously found in SiC>2 
powders at room temperature and comparable at 100 K. 
The high Muonium yield even at low temperatures is an 
important step towards the development of low emittance 
Mu sources for spectroscopy experiments. 
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